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ABSTRACT 

We have shown previously that many of the properties of persistent accretion-powered millisecond 
pulsars can be understood if their X-ray emitting areas are near their spin axes and move as the 
accretion rate and structure of the inner disk vary. Here we show that this "nearly aligned moving 
spot model" may also explain the intermittent accretion-powered pulsations that have been detected in 
three weakly magnetic accreting neutron stars. We show that movement of the emitting area from very 
close to the spin axis to ~ 10° away can increase the fractional rms amplitude from <0.5%, which is 
usually undetectable with current instruments, to a few percent, which is easily detectable. The second 
harmonic of the spin frequency usually would not be detected, in agreement with observations. The 
model produces intermittently detectable oscillations for a range of emitting area sizes and beaming 
patterns, stellar masses and radii, and viewing directions. Intermittent oscillations are more likely in 
stars that are more compact. In addition to explaining the sudden appearance of accretion-powered 
millisecond oscillations in some neutron stars with millisecond spin periods, the model explains why 
accretion-powered millisecond oscillations are relatively rare and predicts that the persistent accretion- 
powered millisecond oscillations of other stars may become undetectable for brief intervals. It suggests 
why millisecond oscillations are frequently detected during the X-ray bursts of some neutron stars 
but not others and suggests mechanisms that could explain the occasional temporal association of 
intermittent accretion-powered oscillations with thermonuclear X-ray bursts. 

Subject headings: pulsars: general - stars: neutron - stars: rotation - X-rays: bursts - X-rays: stars 



1. INTRODUCTION 

Observations made using the Rossi X-ray Timing Ex- 
plorer (RXTE) have led to the discovery of nine neu- 
tron stars in low-mass X-ray binary systems (LMXBs) 
that produce persistent accretion-powered X-ray oscilla- 
tions with periods equal to their millisecond spin peri- 
ods. These accretion-powered millisecond X-ray pulsars 
(APMXPs) are thought to be weakly magnetic stars that 
have been spun up by accr eting angular momentum (see 
lLamb fc Boutloukosl 120081 and references therein). Re- 
cently, three other neutron stars in LMXBs have been 
found to produce accretion-powered millisecond X-ray 
oscillations intermittently. These intermittent APMXPs 
are the focus of this Letter. 

An intermittent accretion-powered X-ray oscilla- 
tion was first dete cted in HETE J1900.1-2455 
(|Gallowav et alj l2007t ). The 377-Hz oscillation disap- 
peared about two months after it was discovered, even 
though the star was still bright. Oscillations at ~ 
376 Hz have r ecently been disco vered in X-ray bursts 
from this star (| Watts et al.ll2009h . A 442-Hz accretion- 
powered X-ray oscillation was serendipitously detected 
in a ~ 500 s interval durin g a longer RXTE observa- 
tion of SAX J1748. 9-2021 ijGavriil et al.ll2007t) . An in- 
dependent analysis of all the currently available RXTE 
data on this X-ray star detected a 442-Hz oscillation 
in 11 intervals during its 2001 and 200 5 outbursts 
(|Altamirano et al.ll2008t iPatruno et al.ll2009f) . The oscil- 
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lation appeared and disappeared on a timescale of several 
hundred seconds. A search of the 1.3 Ms of RXTE data 
available on Aql X-l discovered an accretion-powered X- 
ray oscillation at its 550-Hz spin freq uency for 150 s dur - 
ing the peak of its 1998 outburst (|Caseila et al.ll2008h . 
The spin frequency of Aql X- l was previously k nown 
from its X-ray burst oscillations ijZhang et al.f l998). Due 
to the poor temporal coverage of most LMXBs and the 
sparseness of the intermittent oscillations discovered so 
far, the total number of intermittent APMXPs in the 
Galaxy is unknown. 

The accretion-powered oscillations of the known in- 
termittent APMXPs share several common character- 
istics. All have nearly sinusoidal waveforms and frac- 
tional rms amplitudes < 3%. The second harmonic (first 
overtone) of the spin frequency was rarely detected in 
HETE J1900. 1—2455; when it was, its amplitude was 
<0.4% (jGallowav et all 12007th No second harmonic has 
been detected in the other two in termittent APMXPs , 
with upper l imits ~ 0.4%-0.9% ()Patruno et all 120091 : 
ICasella et al.l [2008) . The amplitudes of the accretion- 
powered oscillations produced by HETE J 1900. 1-2455 
and SAX J1748.9— 2021 seem to increase short l y be- 
fore or after their X-ray bursts dGallowav et alJ 120071 : 
iGavriil et aTll2007l : [AFamirano et al.ll2008D . 

The nearly aligned moving spot model was originally 
proposed to explain the properties of the persistent 
APMXPs (see lLamb et alJl2006l [20071 . l200cl lLamb et all 
120091 hereafter PI). In this model, the X -ray emitting 
areas are close to the spin axis but move around as the 
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Fig. 1. — Fractional rms amplitudes of the full bolometric waveform (left) and its second harmonic component (right) as a function of 
the observer's inclination, for the spot inclinations shown. These waveforms are for two stable, isotropically emitting, 25°-radius antipodal 
spots on the surface of a 1.4Mq star with a radius of 5GM/c 2 , spinning at 400 Hz. 



accretion rate and structure of the inner disk vary. As 
explained in PI, the model also has implications for the 
nuclear-powered oscillations observed in most of these 
pulsars during their thermonuclear X-ray bursts. Here 
we show that the model can also explain the properties 
of the intermittent APMXPs. The reason is that small 
changes in the latitude of the emitting area can cause 
oscillations to appear and disappear if the e mitting area 
is very close to the star's spi n axis (see also lLamb et al.1 
120081 iBoutloukos et all 120081 ). In contrast to PI, which 
presented results for a wide range of spot inclinations, 
this paper focuses on spot inclinations < 15°. 

2. MODELING AND RESULTS 

It is not yet possible to compute the accretion flows 
and X-ray emission of accreting millisecond pulsars from 
first principles, so simplified models must be used, but 
even these models typically have a dozen or more pa- 
rameters. Despite this complexity, it is possible to iden- 
tify trends and determine whether the parameter values 
needed to explain the observations are plausible. This 
is the approach we follow here. Our findings are based 
on the results of several hundred million waveform com- 
putations. This survey greatly extends previous work. 
To summarize these findings concisely, we present exam- 
ples, describe trends, and discuss the parameter ranges 
consistent with the observations. 

To illustrate our findings, we discuss the bolometric 
waveforms produced by is otropic or Comptonized emis - 
sion like that proposed bv lPoutanen fc Gierlihskil (2003) 
from a single circular spot, two antipodal spots, or two 
spots in the same rotational hemisphere on opposite sides 
of the spin axis, as may occur if the s tellar magnetic field 
evolves as described bv IChen et al.1 ()1998f l. These mod- 
els represent the emission from the stellar surface after 
it has been averaged over the time needed to construct 
a pulse profile, which is typically ~ 10 5 -10 6 times the 
~ 1 ms dynamical timescale near the neutron star. 

For conciseness, we feature the representative case of 
spots with radii of 25° on the surface of a 1.4 Mq star 
with a radius of bGM/c 2 (10.3 km) spinning at 400 Hz, 
but we also describe how the results vary when these 
parameters are changed. We assume radiation from the 



stellar surface reaches the observer without interacting 
with any intervening matter. 

We used the ray-tracing code described in PI to calcu- 
late the waveforms produced by emitting spots at various 
rotational latitudes when viewed at different inclinations. 
This code use s the Schwarzschild p lus Doppler (S+D) ap- 
proximation (|Miller fc L amb 1998), which treats the spe- 
cial relativistic Doppler effects (such as aberrations and 
energy shifts) associated with the rotational motion of 
the stellar surface exactly, but treats the star as spher- 
ical and uses the Schwarzschild spacetime to compute 
the general relativistic redshift, trace the propagation of 
light from the stellar surface to the observer, and cal- 
culate light travel-time effects. The code has been val- 
idated by comparing its results with analytical results 
for special cases and with previously reported numerical 
results (see PI). For the stellar models we consider, ro- 
tational di stortion of the star a nd frame dragging can be 
neglected (jCadeau et al.ll2007t) . 

As discussed in PI, pulse amplitudes are continually as 
small as observed only if the emitting areas are near the 
spin axis, so we concentrate on this geometry. Guided by 
the observations cited in Section [IJ we assume that puls- 
ing is detectable if its fractional rms amplitude exceeds 
0.5% and that the second harmonic (first overtone) com- 
ponent of the pulse is detectable if its amplitude exceeds 
0.5%. 

Figure Q] shows some typical features of the oscillations 
produced by spots located within 20° of the spin axis 
(see also Figures 1 and 4 of PI). The fractional ampli- 
tudes are small (< 8% in these examples) and the second 
harmonic is much smaller than the first, except for ob- 
server inclinations near 90°. The amplitude of the first 
harmonic increases linearly with spot inclination for incli- 
nations this small; it also increases linearly with observer 
inclination up to inclinations ~ 20°, but then flattens and 
decreases for observer inclinations > 20°. The amplitude 
of the second harmonic increases monotonically as the 
spot and observer inclinations increase. 

Effect of spot geometry. Insights into the effect of 
the spot geometry on the detectability of oscillations 
can be obtained by comparing the waveforms produced 
by isotropic emission from a single spot, two antipodal 
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spots, and two spots in the same rotational hemisphere 
(spots with the last geometry are expected if neutron 
vortex motion drives both of the star's magnetic p oles 
toward the same rotation pole; see lChen et al.|[T998f) . 

The oscillation produced by a single spot centered 
within 1° of the spin axis is undetectable by observers 
with inclinations < 35° but becomes easily detectable by 
almost all observers if the spot moves > 4° toward the 
spin equator; the second harmonic is undetectable by all 
observers for spot inclinations <8° (see Figure 1 of PI). 

The oscillation produced by antipodal spots is unde- 
tectable for any observer if the spots are within 1° of the 
spin axis but becomes detectable by almost all observers 
if the spots are ~ 12° from the spin axis (see Figure 1). 
For this inclination, almost all observers will see ampli- 
tudes > 2% while observers with inclinations of 20°-70° 
will see amplitudes of 3%-5%. For this spot inclination, 
the second harmonic is < 0.4% and hence undetectable 
for observers with inclinations < 40° (again see Figure 1). 

The oscillation produced by two spots near the same 
rotation pole can also change from being undetectable to 
being easily detectable if the spots move a short distance. 
For example, spots at the same latitude but separated 
in longitude by 160° will produce oscillations that are 
undetectable by all observers if they are within ~ 1° of 
the spin axis but will produce oscillations that are eas- 
ily detectable (amplitudes ~ 2%-4%) by most observers 
(those with inclinations > 60°) if they move ~ 10° further 
from the spin axis. Their oscillations will also be highly 
sinusoidal (second harmonic amplitudes < 1% for all ob- 
servers and less than 0.5% for observers with inclinations 
<65°). 

These results illustrate two important general points. 
First, oscillations can change from being undetectable 
to being easily detectable if the emitting area is near 
the spin axis and moves a small distance away from it. 
Second, emission from two spots near the spin axis can 
produce highly sinusoidal oscillations even if both are 
visible. These conclusions are not significantly affected 
by variations in the spot size and beaming pattern or the 
stellar mass, compactness, and spin rate, for the expected 
ranges of these properties, as we now illustrate (see also 
Figures 2 and 3 of PI). 

Effect of spot size. As discussed in PI, the dependence 
of the fractional modulation on the spot size is relatively 
weak for spots with radii <45°. The reason for this is 
that the width of the radiation pattern produced by spots 
this small is determined primarily by gravitational light 
deflection and the beaming pattern from the stellar sur- 
face. For example, for an observer at an inclination of 
60°, the fractional amplitude of the waveform produced 
by two 5° antipodal spots is only 1.7 times greater than 
the amplitude produced by two 45° spots. The depen- 
dence of the fractional modulation on the spot size is 
even weaker for single spots, increasing by at most 10% 
as the spot radius decreases from 45° to 5°. 

Effect of stellar compactness. As also discussed in PI, 
stars that are more compact tend to produce a lower frac- 
tional modulation, because the stronger light deflection 
averages the flux over a larger fraction of the stellar sur- 
face, reducing the observed flux variation. For example, 
a 1.4M@ star spinning at 400 Hz with two 25° antipodal 
spots 8° from the spin axis will produce a 5% modula- 
tion when viewed at an inclination of 60° if its radius is 



5 GM/c 2 but a 1.5% modulation if its radius is 4 GM/c 2 . 
As another example, a star with the values of the other 
parameters listed above but a radius of 10.5 km will pro- 
duce a 4.2% modulation if its mass is 1.2 M© but only 
a 1.6% modulation if its mass is 1.8 Mq. However, this 
reduction is not in itself sufficient to explain the rare de- 
tection of accretion-powered oscillations in neutron stars 
in LMXBs, because the lowest detectable modulation is 
much lower than these values. Compactness is not a 
plausible explanation of intermittency because it changes 
only very slowly. 

Effect of spin. Stars that are spinning more rapidly 
tend to produce stronger oscillations because their higher 
surface velocities cause larger Doppler boosts and aber- 
rations, making the radiation pattern more asymmetric. 
Fractional amplitudes also tend to increase more steeply 
with spot inclination for stars that are spinning more 
rapidly. The shapes of these curves depend on the spin 
rate and whether the observer sees the spot occulted 
by the star during its rotation. For example, an ob- 
server at 60° sees the fractional modulation from two 
antipodal spots inclined at 8° increase nearly linearly 
from 1.0% at 100 Hz to 4.2% at 600 Hz, whereas an 
observer at 25° sees the fractional modulation from the 
same two spots increase from 2.5% to 3.5% in a non- 
linear fashion, because the spots are occulted for the 
latter observer. The fractional modulation from one or 
two spots increases fairly linearly for spot inclinations 
< 10°, with a slope that increase s with the spin rate. (See 
iPoutanen Sz Beloborodo v 2006 for an analytical approx- 
imation valid when the emitting area is not occulted.) 

Effect of the beaming pattern. The radiation beam- 
ing pattern can affect the amplitude and harmonic con- 
tent of oscillations more than the spot size or stel- 
lar compactness and spin rate, but these effects do 
not change the basic conclusions. To illustrate this, 
we discuss here the waveforms produced by the fan- 
lik e Comptonized emission bea ming pattern proposed 
by IPoutanen &: Gierlihskil (|2003f ). which differ somewhat 
from the waveforms produced by isotropic emission. We 
again consider a single spot, two antipodal spots, and 
two spots in the same rotation hemisphere. 

The oscillation produced by Comptonized emission 
from a single spot is undetectable by all observers with 
inclinations less than 70° if the spot is within 1° of the 
spin axis. If the spot moves to a position 12° from the 
spin axis, the oscillation amplitude increases to 4%, mak- 
ing it easily detectable. 

The modulation produced by Comptonized emission 
from two antipodal spots can be up to twice as large as 
that produced by isotropic emission from the same spots. 
As before, movement of the emitting area in latitude by a 
few degrees will cause the oscillation to change from be- 
ing undetectable to being detectable with an amplitude 
of a few percent, for most viewing directions. 

The oscillation produced by Comptonized emission 
from two spots equidistant from the spin axis and 160° 
apart in azimuth is undetectable by all observers if they 
are both within 2° of the spin axis and is undetectable 
by observers with inclinations less than 45° if they are 
within 8° of the spin axis. Spots with this beaming pat- 
tern can produce modulation fractions as large as the 
largest observed in the APMXPs when they are > 12° 
from the spin axis. 
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The amplitudes and harmonic content of the wave- 
forms produced by these emission geometries and beam- 
ing patterns are similar to those observed in the intermit- 
tent APMXPs. The fractional amplitude of the second 
harmonic is less than 1% for spot inclinations < 12° and 
observer inclinations < 70°, except when there are two 
spots near the same spin pole, in which case the ampli- 
tude of the second harmonic can be comparable to that 
of the fundamental. 

Effect of rapid spot movements. In a more realistic 
model that includes the expected movement of the emit- 
ting area on the surface of the star on timescales as short 
as the ~ 1 ms dynamical time there, the range of incli- 
nations that produce undetectable oscillations is likely 
to be larger than the range discussed above. The rea- 
son is that when the emitting area is near the spin axis, 
its azimuthal motion tends to produce large phase vari- 
ations that, after averaging by the waveform reconstruc- 
tion process, will reduce its apparent amplitude (see PI). 
This effect may help to explain the difficulty of detecting 
accretion-powered oscillations in accreting neutron stars 
with millisecond spin periods. It also means that the 
emitting area may have to move further from the spin 
axis for oscillations to become detectable. 

As the emitting area moves away from the spin axis, 
the emission pattern becomes less axisymmetric and fluc- 
tuations in its azimuthal position produce smaller phase 
variations. Both effects increase the fractional modula- 
tion. 

Motion of the omitting area on timescales shorter than 
the intervals of data used to construct a pulse profile will 
appear as excess background noise (jLamb et all 119851: 
PI). In the physical picture proposed here, the strength 
of this excess noise should be correlated with the ampli- 
tude and harmonic content of the pulse profile (see PI). 
Motion of the emitting area on timescales longer than 
the intervals of data used to construct a pulse profile will 
produce correlated changes in the amplitude and har- 
monic content of the pulse profile that can be measured 
(again see PI). 

3. DISCUSSION 

As explained in PI, the magnetic poles of accreting 
neutron stars are expected to move close to the spin axis 
as the stars spin up to millisecond periods. The resulting 
magnetic field will channel accreting gas toward the axis, 
causing the X-ray emitting area to be close to the spin 
axis. In this Letter we have shown that if the emitting 
area is centered within a few degrees of the spin axis, 
the observed flux modulation at the spin frequency will 
often be so small for a wide range of viewing angles that 
it will be undetectable using current instruments. Rapid 
fluctuations in the position of the emitting area will con- 
tribute to making the oscillation at the spin frequency 
undetectable. These effects may explain why accretion- 
powered oscillations have have been detected in so few 
accreting neutron stars with millisecond spin periods (see 
also PI). 

We have shown further that if the emitting area wan- 
ders toward the spin equator by ~ 10°, the oscillation 
can change from being undetectable to being easily de- 
tectable, with amplitudes comparable to those observed 
in the intermittent APMXPs. In most cases, the am- 
plitudes of the second and higher harmonics would be 



detectable only when the amplitude of the fundamental 
exceeds several percent. 

Our computations show that this behavior occurs for 
a variety of spot geometries and radiation beaming pat- 
terns, and a wide range of spot radii, stellar masses and 
radii, and observer inclinations. Infrequent episodes of 
detectable accretion-powered oscillations at the spin fre- 
quency are more likely for stars that are more compact, 
because such oscillations will generally have lower am- 
plitudes in these stars and are therefore less likely to be 
detectable. 

In addition to providing a possible explanation for the 
sudden appearance of accretion-powered oscillations at 
the spin frequency, our model predicts that otherwise 
persistent oscillations may occasionally disappear if the 
center of the emitting area moves too close to the spin 
axis. Such oscillation dropouts, which could occur on a 
variety of timescales, have not yet been reported, but 
may be present in the existing data on known APMXPs 
or found in data on APMXPs discovered in the future. 
We encourage searches for such dropouts. 

The basic physical picture of accreting millisecond X- 
ray pulsars proposed here and in PI may also explain the 
observed pattern of millisecond oscillation detections in 
the X-ray bursts of different types of neutron stars. 

Assume for the sake of argument that the mag- 
netic poles of most weak-field accreting neutron stars 
in LMXBs are so closely aligned with their spin axes 
and the accretion flow so symmetric around the rotation 
pole that accretion-powered X-ray oscillations are unde- 
tectable (see PI). If the surface magnetic field is strong 
enough to create a symmetric pattern of nuclear burning 
around the spin axis (again see PI), X-ray burst oscilla- 
tions will also be undetectable. 

In PI we conjectured that the magnetic fields 
of the persistent APMXPs SAX J1808.4-3658 and 
XTE J1814— 338 are slightly asymmetric and strong 
enough to confine nuclear fuel, causing their nuclear 
burst oscillations to be locked or nearly locked to their 
accretion-powered oscillations, whereas the magnetic 
fields of other persistent APMXPs are too weak to force 
this. 

The model of the intermittent APMXPs proposed here 
also suggests why accretion-powered oscillations and X- 
ray bursts may be associated in ti me, as appears to 
be the case in SAX J1 748. 9-2021 (|Altamirano et all 
2008; Pa truno et aHl2009f ). The appearance of accretion- 
powered oscillations signals a shift in the impact point of 
the accreting matter. On the one hand, such a shift could 
be due to a change in the accretion flow through the inner 
disk. This is likely to change the thermal structure and 
fuel loading in the outer layers of the neutron star, which 
might trigger a nuclear burst. On the other hand, a nu- 
clear burst is expected to temporaril y change the accre- 
tion flow in the inner disk (see, e.g.. lMiller fc Lamb| [i996 
and references therein) , changing the impact point of the 
accreting matter and thereby possibly causing accretion- 
powered oscillations to become detectable. 
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